EXPERIMENTAL  AVAILABILITY  TABLES 
FOR  FINITE  SPARES  BACKLOGS 


KH  Ju  Park 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 

EXPERIMENTAL  AVAILABILITY  TABLES 
FOR  FINITE  SPARES  BACKLOGS 

by 

Park,  Kil  Ju 

March  1979 

Thesis 

Advisor:                 J.  D. 

Esary  ", 

Approved  for  public  release;  distribution  unlimited 


T188632 


security  classification  of  Tmi  p»ct  fjWwg  d«»  Entarad)  phouyknoxu. 


WaTBffl 


REPORT  DOCUMENTATION  PAGE  ««w«^5miS> 


I.     REPORT   KUMtt* 


•WOttL  READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 


2.   OOVT    ACCESSION   MO, 


1.     RECIPIENT'S  CAT  ALOG   NUMBER 


4.     TITLE   (and  Subilllm) 

Experimental  Availability  Tables 
for  Finite  Spares  Backlog 


S.     TYPE  OF    REPORT   a   PCRIOO  COVERED 

Master's  Thesis; 


March  1979 


•  ■  PERFORMING  ORG.  REPORT  NUMICR 


7.     AUTHORfaJ 


t.     CONTRACT  OR  GRANT  NUMSERfaj 


Park,   Kil  Ju 


t.     PERFORMING  ORGANIZATION   NAME  AND  AOORESS 

Naval  Postgraduate  School 
Monterey,  California   93940 


10.     PROGRAM   CLEMENT,  PROJECT,    TASK 
AREA   •   WORK   UNIT  NUMBERS 


It.     CONTROLLING  OFFICE  NAME   ANO   AOORESS 

Naval  Postgraduate  School 
Monterey,  California   93940 


12.     REPORT  DATE 

March   19  79 


IS.     NUMIIK  OF   PAGES 


36 


14.     MONITORING   AGENCY   NAME  *    AOORESV"  dlltatwnt  ham   Control Unv   Olll cm) 


IS.     SECURITY   CLASS,  (at  thla  riport) 

Unclassified 


11a.    DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


l«.     DISTRIBUTION  STATEMENT  (ol  thla  Hamart) 


Approved  for  public  release;  distribution  unlimited 


17.     DISTRIBUTION  STATEMENT  (al  lha  moan-act  mntatad  In  Black  30,  II  dlUatant  /ran  Raparrt) 


It.     SUPPLEMENTARY  MOTES 


I*.     KEY  WORDS  (Contlmia  on  tawawaa  ml  am  II  nmcaaamry  anal  Identity  ay  aloe*  nnavj 

operational  availability 
experimental  tables 


20.     ABSTRACT  (Conittwa  an  rawataa  atda  II  nacaaaary  and  Idanttty  ay  alack  mamavt) 

Experimental  tables  of  availabilities  at  time  t  are  obtained 
for  a  device  whose  performance  is  described  by  an  alternating 
renewal  process  with  a  finite  number  of  failure-renewal  cycles, 
corresponding  to  having  a  finite  spares  backlog.   Failure  and 
repair  rates  are  assumed  to  be  constant,  and  attention  is 
restricted  to  cases  in  which  the  repair  rate  is  larger  than 
the  failure  rate. 


DD    fom" 


I  JAN  71   1473 

(Page   1) 


EDITION  OF   I  NOV  •■  IS  OBSOLETE 

S/N    0103-01*-  *601    | 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Data  gntarad) 


Approved  for  public  release;  distribution  unlimited 

Experimental  Availability  Tables 
for  Finite  Spares  Backlogs 


by 


Park,  Kil  Ju 

Lieutenant /'Republic  of  Korea  Navy 

B.S.,  Republic  of  Korea  Military  Academy,  1973 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  OPERATIONS  RESEARCH 


from  the 

NAVAL  POSTGRADUATE  SCHOOL 
March  19  79 


i»5 


ABSTRACT 

Experimental  tables  of  availabilities  at  time  t  are 
obtained  for  a  device  whose  performance  is  described  by  an 
alternating  renewal  process  with  a  finite  number  of  failure- 
renewal  cycles,  corresponding  to  having  a  finite  spares 
backlog.   Failure  and  repair  rates  are  assumed  to  be  constant, 
and  attention  is  restricted  to  cases  in  which  the  repair  rate 
is  larger  than  the  failure  rate. 
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I.   INTRODUCTION 

The  most  commonly  encountered  working  definition  of  the 
"availability"  of  a  device, 

MTTF 


(1.1)  Availability  = 


MTTF  +  MTTR 


represents  the  long-term  or  steady-state  probability  that 
the  device  will  be  found  in  an  "up"  or  functioning  condition 
when  two  specific  conditions  are  satisfied.   One  condition  is 
that  there  is  an  alternation  of  failure  and  repair  cycles 
in  which  times  to  failure  and  times  to  repair  are  independent 
realizations  from  some  failure  and  repair  distributions 
satisfying  minimal  regularity  conditions.   The  second,  and 
here  most  important  condition,  is  that  the  alternation  of 
failure  and  repair  continues  indefinitely,  so  that  the 
performance  of  the  device  is  described  by  a  standard  alternating 
renewal  process. 

For  many  equipments,  the  second  condition  cited  above 
implies  access  to  an  infinite  backlog  of  spares.   In  many 
operational  contexts  this  sort  of  spares  support  cannot  be 
realistically  assumed. 

If  availability  is  considered  to  derive  from  an  alter- 
nating renewal  process  with  a  finite  number  of  cycles, 
corresponding  to  a  finite  backlog  of  spares,  then  expressions 
for  availability  become  complex  as  compared  with  equation  (1.1). 


This  thesis  is  devoted  to  computational  experiments  with 
some  "finite  spares"  availability  expressions.   The  end 
objective  of  such  experiments  is  to  be  able  to  determine  the 
circumstances  in  which  equation  (1.1)  furnishes  an  adequate 
approximation,  or  alternately  to  be  able  to  provide  compu- 
tationally feasible  alternatives  to  its  use. 


II.   MATHEMATICAL  MODEL 

The  mathematical  model  on  which  the  usual  expressions 
for  the  availability  of  a  device  are  based  is  an  alternating 
renewal  process  with  an  infinite  number  of  failure-repair 
cycles.   In  situations  where  repair  requires  replacement  of 
a  failured  device  by  a  spare,  this  corresponds  to  having  an 
infinite  number  of  spares.   The  model  studied  here  is  modi- 
fied to  allow  only  a  finite  number  of  failure-repair  cycles, 
corresponding  to  having  a  finite  number  of  spares. 

The  simplest  assumptions  about  failure  and  repair  times 
are  made;  failure  rates  are  constant,  and  repair  rates  are 
constant.   Only  those  processes  that  begin  with  a  functioning 
device  installed  are  considered. 

In  greater  detail,  the  failure-repair  process  considered 
is  as  shown  in  figure  2.1, 


n 


end 


Figure  2.1   Failure-repair  process 


where  n  is  the  number  of  spares ,  Tfl  is  the  time  to  failure 
for  the  original  device,  T,,  T  , T   are  the  times  to 

failure  for  the  n  spares,  and  Sn,  S. , S   ,  are  the 

r    '      0 '   1 '       n-1 

times  to  replace  the  original  device  and  the  first  n-1  spares 

It  is  assumed  that  Tn ,  Sn,  T   , ,  S   , ,  T   are  inde- 

0 '   0 '        n-1 '   n-1 '   n 

pendent  random  variables,  and  that  Tfl, ,  T  are 

exponentially  distributed  with  failure  rate  X,   while 

S„, ^n-1  are  exPonent:ially  distributed  with  repair 

rate  r\ . 

The  availability  at  time  t,  A  (t)  ,  of  the  original 
device,  supported  by  its  backlog  of  n  spares,  is  the 
probability  that  the  process  shown  in  Figure  2.1  is  in  an 
"up"  condition  at  time  t;  i.e.,  that  at  time  t  either  the 
original  device  or  one  of  its  spares  is  installed  and  still 
functioning. 

The  increment  in  availability  at  time  t  due  to  the 
k   spare,  I,  Ct)  is  defined  by 

(2.1)       Ik(t)  «  Ak(t)  -  Ak_1(t)    k  =  1, ,  n. 

Before  proceeding  to  a  derivation  of  expressions  for 

I  (t)  and  A  (t)  in  a  general  case,  two  special  cases  are 
n        n        °  '  r 

considered;  repair  rate  n  equal  to  infinity,  and  repair 
rate  n  equal  to  failure  rate  X.   These  are  boundary  cases 
for  the  cases  of  likely  practical  interest,  in  which  it  is 
reasonable  to  expect  that  repair  rate  will  exceed  failure 
rate. 


In  any  case,  AQ(t) ,  availability  at  time  t  with  no 
spares  is  given  by 

(2.2)         AQ(t)  =  P[TQ  >  t]  =  e"At,  t  >  0. 

In  the  following  sections,  it  will  be  convenient  to  let 
(2.3) 


u     = 
n 

To 

+ 

,  .    +   T    , 
n' 

Vn   - 

so 

+ 

n' 

w     = 
n 

u 

n 

+ 

V 
n 

=  CS0  +  V  +  +  <Sn  +  V 


A.   REPAIR  RATE  EQUAL  TO  INFINITY 

The  simplest  case  is  the  one  in  which  no  time  is  required 
to  repair  a  failed  unit,  provided  a  spare  unit  is  available. 

In  this  case  the  contribution  of  the  first  spare  is 


(2.4)       I  Ct)  =   P[T   >  t-s|UQ  =  s]  fy  (s)ds 

l0  l         o 


where 


f   (s)  =  Xe'Xs ,  s  >  0 
u0 


is  the  gamma  {1,X}  density.   Thus 

(2.5)  Ix(t)  =  I   e"X(t"s)Ae"Xs  ds 

Jo 


,  -As 
=  Ae    t, 


and  the  availability  of  a  system  having  one  spare  is 


(2.6)  A^t)  =  AQ(t)  +  I^t) 

=  e"Xt(l  +  At) 


The  contribution  of  the  second  spare  is 


C2.7)      I2(t)  =  |   P[T2  >  t-s|U1  =  s]fu    (s)  ds  , 

0  1 


where 


,2      -As 

Vs5  -  t-ftzr  's-°' 


is  the  gamma  {2, A}  density.   Thus 


(2.8)  I2(t)  =  |   e-X^-s^  xZf,^j5  ds 


,,«..2  -At  1 

=  (At)  e    2T 


and  the  availability  of  a  system  having  two  spares  is 

2 
(2.9)  A2(t)  =  e"Xt(l  +  At  +  y$—)     • 

Generally,  the  contribution  of  the  n   spare  is 


10 


(2.10)         I    ft)    =    (      P[Tn   >    t-s|Un_1   =    sjfy        (s)    ds    , 

Jo.  n_1 


where 


,n  n-1    -As 

£u    .  C«)  -      srUj        >  5  >  °  > 

n-l  v   J 


is   the   gamma    {n,A}   density.      Thus 


(2.11)  I.,,  Ct)    -  e  As        e  , 

n  i  ^ as 


-    (Xt)      e  — ,    , 


and  the  availability  of  a  system  having  n  spares  is 
(2.12)     An(t)  =  e"Xt(l  +  At  +  I££L  +  +  ^   ) 


Availabilities  A  (t)  obtained  from  equation  (2.12)  are 
shown  in  table  1.   In  this  table  n  represents  the  number  of 
spares,  and  the  contribution  I  (t)  of  the  n   spare  can  be 
found  by  subtracting  A  _-<(t)  from  A  (t) .   The  availabilities 
shown  in  the  last  column  are  for  an  alternating  renewal  pro- 
cess with  an  infinite  number  of  failure-repair  cycles.   This 
corresponds  to  having  an  infinite  number  of  spares. 

B.   REPAIR  RATE  EQUAL  TO  FAILURE  RATE 

The  failure-repair  process  considered  is  that  in  which 
failure  times  and  repair  times  are  exponentially  distributed 
with  equal  rates,  i.e.,  X   =   ri- 
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In  this  case,  the  contribution  of  the  first  spare  is 


rt 


(2.13)  I^t)    =        P[TX    >    t-s|WQ    =   s]fw    (s)    ds    , 

'o  ° 


where 


Vs)  =  TTO-"  's-0' 


is  the  gamma  {2, A}  density.   Thus 


(2.14)         Ix(t)  =  J   e-A(t-s3  \2^J"XS  ds 

0 


, 2  -At  .2 
a  e    t 

"  r(2)  T"  ' 

and  the  availability  of  a  system  having  one  spare  is 


Ax(t)  =  AQ(t)  +  IL(t) 


=  e^  (1  «■  ^    . 


Generally,  the  contribution  of  the  n   spare  is 


(2.15)    In(t)  =  (   P[Tn  >  t-s|Wn_1  =  s]fw    (s)  ds  , 

0  n~1 


where 


->  2n  2n-l  -As 

£W   ,  &    =     kin)"  •  S  i  °  ' 

n-1 


13 


is    the   gamma   {2n,A}    density.      Thus 


(2.16) 


I    (t)    = 

nv    J 


A(t-s)    Xe 


As 


TTTn) 


,-,  „ -.  2n- 1    j 
(As)  ds 


X2ne-Xt    f      2n-l    , 
(2n-lj!  S  ds 

J  0 


(Xt)2ne-Xt 
"         (2nj! 


and   the   availability  of   a   system  having  n  spares    is 


(2.17)        An(t)    =   e"At    (1    +  l£fL.  + 


C.   REPAIR  RATE  GREATER  THAN  FAILURE  RATE 

The  failure-repair  process  considered  in  this  section 
is  visualized  as  one  in  which  the  repair  rate  is  greater 
than  the  failure  rate.   This  influences  the  format  in  which 
the  results  are  displayed. 

■f"Vi 

In  this  case,  the  contribution  of  the  n   spare  is 

n 


(2.18)   In(t)  = 


P[T„  >  t-s|Wn_1  =  s]  f     (s)  ds  , 


0 


n-1 


where 


£W     CS)  -  /    fy 

n-1     J  0   n-1 


(s)  £     (t-s)  ds  , 


■u 


n-1 
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and 


n-1  ^   J 


is    the   gamma    {n,n}   density, 
while 


n 


fn        Ct-s)    =  jrr-r  (t-s) 


■u 


n-1 


TJnJ 


n-1      -A(t-s) 


,    t-s     •    0    , 


is    the   gamma    {n,A}   density.      Thus 


(2.19)    In(t)    =    |      e 


t  C   S  r, 

A Ct-s) 


J 


n 


0 


n ..n-l-nu'     X r_   .^n-l 

rTnT  U       e  TJnJ   (s_u) 


-X(t-S)       ,        j 

•    e      v        ;    du   ds 


Inverting  the  order  of  integration,  equation  (2.19) 
becomes 

<t  nt 


(2.20)    In(t)    = 


0  ^u 


-X  (t-s)      ri11     ,,n-l    -nu      Xn      ,        ^n-1 
e  TTHTU       e  TJnJ  Cs"u) 


-X(s-u)     ,e 
*    e      *•        J    ds 


du 


_      -xt       nnxn      I        n-1   -(n-X)u 
=    e  — — «■    I      u       e 


rt 


{r(n)}"  J0 


(s-u)  ds   du 


J  u 


Let   v  =   s    -   u.      Then 


0 


(s-u)11"1   ds    = 


r 


t-u 


vF-1   dv  =    (t-u) 


n 


n 
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Thus  equation  (2.20)  reduces  to 


(2.21)  In(t) 


=  e 


-xt  nV1 
rcn)r(n+u 


.,n-l ,.   .n  -  (n-X)  u  , 
u   (t-u)  e  v   J      du  , 


and  the  availability  of  a  system  having  n  spares  is 
(2.22)      An(t)  =  AQ(t)  +  Ix(t)  +  +  In(t)  . 

Note  than  when  r\    =    X ,    equation  (2.21)  reduces  to 
equation  (2. '6),  since  then 


ln& 


=  e 


•Xt     A 


2n 


i'ujru+ij 


^^(t-u)11  du 


Let  u :.«■  tv.   Then 


r 


un_1(t-u)n  du 


r 


(tv)n"1{t(l-v)}n  tdv 


.  2n    n-1 ,,   ->  n  , 
t     v    (1-v)   dv 


_  ^2n  r(n)T(n+l) 


=  t 


T(2n+1) 
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so    that 


I    (t) 


\t    2n 
e      z\  2n   r(n)r(n+l) 

r(njr(n+lj    Z  i'C2n+lJ 


(At)2VXt 
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III.   APPROXIMATION  OF  THE  MATHEMATICAL  MODEL 

In  section  II,  we  derived  a  mathematical  model  for 
availability. 

In  this  section,  we  discuss  methods  of  approximation  to 
obtain  numerical  values  of  availability. 

A.   EXPONENTIAL  EXPANSION  APPROXIMATION 

The  integral  in  equation  (2.21)  can  be  approximated  by 
expanding  its  exponential  term,  i.e., 


(3.1) 


Cn-X)   ,   .  n,  4  (n-x)   2 
k    =  1  -  (n-A.)u  +  ^  '  ■  i —  u 


2T 


Thus  the  integral  becomes 


r 


(3.2) 


\P'ltt-u)n   {i-(n-x)u+Iil^l- ^2 


}  du 


un-1(t-u)n  du 


r 


-   (n-X) 


u^(t-u)n  du 


2  r* 


u    (t-u)   du 


We  know  that , 
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u"_J-(t-u)"  du  =  t 


.   „n-l^_„,n  An    m    «.2n  r(n)r(n+l) 


r(2n+l) 


Thus  equation  (3.2)  becomes 

f3  31       I  ftl  «  e"Atn^n  ,  ft2n  r(n)r(n+l) 

C3-3J       VtJ    r(n)r(n+l)  U    r(2n+l) 

rn  xit2n+1  r(n+pr(n+i) 


+    (n-A)2t2n+2  r(n+2)r(n+l) 
21 l'(2n+3J 


e  ■ 


(nt)"(xt)"  ,   r(n)  ,   .    Xl0   r(n+i) 

rCnJ 1r(2n+lJ    "    lnt   AtJ    r(2n+2) 


.    (nt-xt)2  r(n+2) 
2~!  f  (,2n+  3  J 


} 


Computational  experiments,  with  the  approximation 
represented  by  equation  (3.3)  have  indicated  unsatisfactory- 
convergence  behavior  when  (nt-Xt)  is  large,  a  case  of  some 
practical  interest,  and  so  this  approach  was  not  pursued. 

B.   SIMPSON'S  RULE  APPROXIMATION 

The  integral  in  equation  (2.21)  can  be  approximated  by 
using  Simpson's  rule. 


20 


Let  v  =  £  •   Then 


(3.4) 


u    (t-u)e  *•   ;u  du 


n-1 


CtV)    -{t(l-v)}n  e-(n'A)tVt 


dv 


=  t2n 


r1 


V 


n"1Ci-v)n  e-Cnt-\t)v  dy 


Thus  equation  (2.21)  becomes 

-xt  (nt)nqt)n  f 

(3.5)      In(t)    =   e    At   I'CnJl'tn+lJ  J 


n,,   ,n   -  (nt-Xt)  v 
v  (1-v)   ©         dv 


Now  Simpson's  rule  can  be  applied  to  the  integral  in 
equation  (3.5)  to  obtain  numerical  values  of  availability. 
Simpson's  rule  as  applied  is 


r 


(3.6) 


J 


f(x)dx  =  3(y1+4y2H-2y3  +  4y4+....+2yni_1  +  4ym+ym+1) 


where 


h  = 


b-a 
m 


and 


X1  =  a,  X2  =  a+h,  ,  Xm+1  =  a+mh  =  b, 
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while  y1   =    f(X1),    Y2   =    fC*2)  »    '    ym+l   "    ^m+l5 
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IV.   TABLE 

In  this  section  experimental  tables  of  availabilities 

A  (t)  are  shown,  which  were  obtained  from  the  mathematical 
n 

model  evaluated  using  Simpson's  rule,  i.e.,  using  equation  (3.6) 
to  evaluate  equation  (3.5)  with  m  =  500,  X,  =  0.0001  and 
^"501  =0-9999.   The  reason  for  choosing  m  =  500  is  that  compu- 
tational experiments  with  choices  of  m  greater  than  or  equal 
to  500  gave  a  stable  and  accurate  result. 

In  these  tables  n  represents  the  number  of  spares, 
and  the  contribution  of  the  n   spare  I  (t)  can  be  found  by 
subtracting  A  _,(t)  from  A  (t) . 

The  availabilities  shown  in  this  section,  Table  3-12, 
are  for  cases  in  which  nt  >  At. 

The  availabilities  shown  in  the  last  column  are  for  an 
alternating  renewal  process  with  an  infinite  number  of  failure- 
repair  cycles .   This  corresponds  to  having  an  infinite  number 
of  spares. 
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V.   SUMMARY  AND  CONCLUSIONS 

Certain  computational  approaches  were  tried  for  obtain- 
ing availabilities  for  a  device  supported  by  only  a  finite 
backlog  of  spares,  using  the  simple  assumptions  that  failure 
and  repair  rates  are  constant.   Real  failure  and  repair  distri- 
butions may  be  more  complex,  but  the  case  considered  is  a 
good  case  for  initial  computational  experiments. 

Of  the  two  approaches  tried,  neither  proved  entirely 
satisfactory  in  obtaining  availabilities  in  a  way  that  is 
fast  and  suitable  for  use  with  small-scale  computational 
facilities,  e.g..  hand-held  calculators.   Also  neither  was 
effective,  over  the  entire  range  of  failure  rate  and  repair 
rate  combinations  that  might  be  of  interest. 

Since  an  easily  used,  readily  accessible, way  to  assess 
the  impact  of  finite  spares  backlogs  on  availability  is 
desirable  in  many  mission  planning  contexts,  further  computa- 
tional approaches  should  be  tried. 

The  tables  presented  in  section  IV  give  availability 
values  with  which  the  result  of  such  experiments  can  be 
compared. 
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